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19.  ABSTRACT  (continued) 


A  two-dimensional  computer  model  was  developed  to  predict  the  thermally  induced  stress 
in  a  growing  Czochralski  GaAs  boule  from  fundamental  growth  parameters.  A  "stress-conversion 
coefficient4'  was  obtained  from  experimental  data  and  was  applied  to  the  stress  results  to 
predict  two-dimensional  EL2  profiles.  These  predicted  profiles  were  compared  with  experi¬ 
mental  results, 

A  quasi-two-dimensional  annealing  computer  model  was  constructed  to  model  diffusion  from 
calculated  and  experimental  EL2  distributions.  Experimental  annealing  data  was  used  in  the 
annealing  model  to  obtain  the  approximate  diffusion  constant  ol  the  lil.?  delect.  This 
diffusion  constant  was  then  used  in  conjunction  wllh  the  annealing  model  to  predict  the 
theoretical  LL2  concentration  before  and  after  annealing.  The  predicted  profiles  were 
compared  with  experiment. 

A  Fermi-level  computer  model  was  constructed  to  determine  the  Fermi-level  at  each  point  in 
the  wafer  as  a  function  of  four  arbitrary  impurities.  Once  obtained,  the  Fermi  level  was  used 
to  calculate  the  resistivity  at  each  location  in  the  wafer.  Theoretical  results  were  calcula¬ 
ted  for  unannealed  material  containing  EL2,  carbon,  iron  and  silicon. 
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I.  REVIEW 


The  primary  objective  of  this  process  modeling  program  was 
to  identify  Semi-Insulating  Gallium  Arsenide  (SI-GaAs)  substrate 
and  implantation  parameters  that  affect  device  behavior. 

Three  major  areas  were  investigated. 

1.  Predicting  the  two-dimensional  distribution  of  the 
EL2  anti-site  defect. 

2.  Modeling  changes  in  the  two-dimensional  EL2 
distribution  resulting  from  the  annealing  of  SI-GaAs. 

3 .  Determining  the  electrical  characteristics  of  the 
SI-GaAs  substrate  and  the  thermal  stability  of  these 
characteristics . 

A  two-dimensional  computer  model  was  developed  to  predict 
the  thermally  induced  stress  in  a  growing  Czochralzki  GaAs  boule 
from  fundamental  growth  parameters.  A  "stress-conversion 
coefficient"  was  obtained  from  experimental  data  and  was  applied 
to  the  stress  results  to  predict  two-dimensional  EL2  profiles. 
These  predicted  profiles  were  compared  with  experimental  results. 


A  quasi-two-dimensional  annealing  computer  model  was 
constructed  to  model  diffusion  from  calculated  and  experimental 
EL2  distributions.  Experimental  annealing  data  was  used  in  the 
annealing  model  to  obtain  the  approximate  diffusion  constant  of 
the  EL2  defect.  This  diffusion  constant  was  then  used  in 
conjunction  with  the  annealing  model  to  predict  the  theoretical 
EL2  concentration  before  and  after  annealing.  The  predicted 
profiles  were  compared  with  experiment. 

A  Fermi-level  computer  model  was  constructed  to  determine 
the  Fermi-level  at  each  point  in  the  wafer  as  a  function  of  four 
arbitrary  impurities.  Once  obtained,  the  Fermi  level  was  used 
to  calculate  the  resistivity  at  each  location  in  the  wafer. 
Theoretical  results  were  calculated  for  unannealed  material 
containing  EL2,  carbon,  iron  and  silicon. 
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II.  BACKGROUND  ON  SI-LEC  GaAs  MATERIAL 

1.  Introduction 

Semi-insulating  Gallium  Arsenide  (SI-GaAs)  substrates  offer 
the  advantages  of  device  isolation,  low  interconnect  capacitance, 
simplified  processing,  and  direct  ion- implantation  for  device 
fabrication.  However,  the  quality  and  homogeneity  of  SI-GaAs 
material  varies  widely  due  to  difficulties  encountered  during 
crystal  growth.  In  order  to  achieve  the  tremendous  advantages 
of  GaAs  technology,  it  is  essential  to  understand  the  physical 
mechanisms  behind  the  growth  of  the  GaAs  material. 


2 .  Growth  methods 


Two  major  growth  techniques  (Horizontal  Bridgeman  and 
Liquid-Encapsulated-Czochralski)  exist  for  producing  SI-GaAs 
material.  In  the  Horizontal  Bridgeman  (HB)  method,  the  crystal 
is  grown  by  zone-refining  a  GaAs  melt  in  a  quartz  ampoule.  The 
use  of  quartz  introduces  silicon  as  an  impurity,  requiring 
electrical  compensation  with  chromium  in  order  to  maintain  semi- 


insulating  properties, 


Unfortunately,  chromium  tends  to 


outdiffuse  during  thermal  processing,  causing  the  surface  to 
change  character  from  n  to  p  type  during  thermal  annealing.  This 
thermally  unstable  behavior  means  that  HB  crystals  are 
unacceptable  for  any  annealed  ion-implantation  process  (1]. 

In  the  Liquid-encapsulated-Czochralski  (LEC)  method,  the 
crystal  is  grown  by  a  Czochralski  process  incorporating  an 
encapsulating  blanket  of  boric  oxide  to  maintain  proper 
stoichometry.  One  major  advantage  of  the  SI-LEC  growth  process 
is  that  the  material  can  be  grown  semi-insulating  without 
intentional  doping.  The  unintentional  compensation  mechanism 
responsible  for  the  semi-insulating  behavior  of  LEC-SI  GaAs  is 
believed  to  be  EL2  deep  donors  compensating  carbon  shallow 
acceptors  [2]. 
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The  major  disadvantage  of  the  SI-LEC  process  is  the  high 
dislocation  density  obtained  during  growth  (ten  to  fifty  times 
that  of  the  HB  process)  [3].  The  origin  of  these  dislocations  is 
not  fully  understood,  but  it  appears  to  be  related  to  abrupt 
temperature  gradients  occurring  during  growth.  However,  s  e 
results  on  isoelectronic  doping  of  SI-LEC  material  with  gallium 
and  indium  show  that  impurity  hardening  by  isoelectronic  doping 
can  greatly  reduce  the  dislocation  density  without  significantly 
changing  the  electrical  properties  [4,5] 
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3.  The  Unintentional  Compensation  of  EL2  with  Carbon 


A  crucial  element  in  the  growth  of  SI-LEC  GaAs  is  the 
unintentional  compensation  of  the  EL2  defect  by  carbon 
impurities . 

Carbon  is  probably  introduced  into  the  melt  by  the  graphite 
heaters  used  to  heat  the  crucible.  Carbon  (with  a  segregation 
coefficient  of  approximately  2)  is  constant  across  each  wafer 
slice,  but  decreases  from  the  seed  to  the  tail  of  the  boule  [6]. 

EL2  (a  deep  level  trap  associated  with  the  As/Ga  anti-site 
defect)  increases  from  seed  to  tail  in  an  As-rich  boule, 
decreases  from  seed  to  tail  in  a  Ga-rich  boule  and  is  essentially 
constant  in  a  stoichiometric  boule  [7].  However,  since  EL2  is 
associated  with  As  on  a  Ga  site,  it  is  present  only  in  very  low 
concentrations  in  Ga-rich  boules. 

The  major  concern  with  the  EL2  defect  is  its  variation  in 
the  radial  direction.  As  can  be  seen  in  Figure  1,  the  EL2 
concentration  displays  a  characteristic  "W"  shape  when  measured 
across  a  diameter  [ '?  :  .  This  in  of  great  concern,  because  the 
compensation  between  a  constant  carbon  concentration  and  a  "W"- 
shaped  EL2  concentration  results  in  non-uniformity  of  electrical 
parameters  across 


a  v/afer  [9], 


[EL2J  x  1015cm-3 


16  _ 


14  _ 


12  _ 


rho  =  0 


Normalized  radius,  rho  =  r/r0 


Figure  1  -  EL2  Variation  in  the  Radial  Direction 


The  major  concern  with  the  EL2  defect  is  its  variation  in  the 
radial  direction.  As  can  be  seen  above,  the  EL2  concentration 
displays  a  characteristic  "W”  shape  when  measured  across  a 
diameter.  This  is  of  great  concern,  because  the  compensation 
between  a  constant  carbon  concentration  and  a  "W'-shaped  EL2 
concentration  results  in  non-uniformity  of  electrical  parameters 
across  the  wafer. 
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4 .  Sunnary 


Clearly,  the  key  technical  issue  in  the  growth  of  SI-LEC 
GaAs  is  the  nature  of  the  EL2  defect.  Answers  to  the  following 
questions  are  crucial  to  understanding  the  EL2  defect  and  its 
impact  on  device  electrical  performance.  The  remainder  of  this 
report  will  address  some  of  these  questions. 

1.  What  is  the  physical  mechanism  which  creates  EL2? 

2.  Can  the  EL2  concentration  profile  be  predicted  from 
fundamental  information? 

3 .  What  is  the  origin  of  the  ubiquitous  "W"  shape  for  radial  EL2 
density  profiles? 

4.  Can  the  "W"  shaped  EL2  concentration  profile  be  made 
constant? 

5.  Does  the  compensation  between  the  "W"  shaped  EL2 
concentration  and  the  constant  carbon  concentration  result  in 
variation  in  electrical  parameters? 

6.  Is  the  compensation  between  EL2  and  carbon  thermally  stable? 
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III.  DETERMINING  THE  TWO-DIMENSIONAL  DISTRIBUTION  OF  THE  EL2 
ANTI-SITE  DEFECT. 

1.  Background  -  Identification  and  Evidence  for  the  EL2 
defect 

Little  was  understood  about  the  material  growth  of  GaAs 
until  1981,  when  the  materials  groups  at  both  Rockwell  and 
Westinghouse  began  to  seriously  address  the  material  problem. 

These  researchers  discovered  that  extremely  high  quality 
undoped  GaAs  crystals  could  only  be  reproducibly  grown  from  LEC 
material  if  the  melt  were  maintained  with  an  arsenic  atomic 
fraction  between  0.47  and  0.52.  They  also  observed  that  high 
Hall  mobility  could  only  be  maintained  if  the  melt  were  slightly 
arsenic  overrich.  Finally,  they  noted  that  the  tails  of  gallium 
rich  boules  underwent  a  transition  from  SI  to  p-type  down  the 
length  of  the  boule,  while  arsenic  rich  boules  did  not  display 
this  property  [10]. 

The  conclusion  was  that  a  balance  had  been  reached  between 
shallow  acceptors  due  to  impurities  and  some  sort  of  deep  donor 
which  was  related  to  the  arsenic  concentration.  A  deep  donor 
with  this  type  of  concentration  dependent  behavior  is  likely  to 
be  a  point  defect.  If  one  assumes  that  the  defect  is  a  simple 
atomic  disorder,  then  only  an  anti-site  defect  can  account  for 
the  swing  from  p  to  n  as  the  material  becomes  arsenic  rich  [11’ . 

9 


A  big  step  in  identifying  the  principle  defect  in  GaAs 
material  was  taken  by  Martin  et.  al.  In  1977,  he  had  summarized 
all  of  the  currently  known  defects  in  GaAs  as  measured  by  DLTS 
(Deep  Level  Transient  Spectroscopy)  as  well  as  their  cross- 
sections  and  activation  energies  [12].  In  1980,  Martin  et.al. 
used  OTCS  (Optical  Transient  Current  Spectroscopy)  to  identify 
the  deep  donor  they  named  EL2,  with  an  activation  energy  of  about 
0.784  eV,  and  was  identified  with  a  characteristic  optical 
absorption  fingerprint  between  1  and  1.4  microns.  [13] 

At  this  point  in  time  (1981),  the  conclusion  was  that  an 
anti-site  defect  existed  in  the  high-quality  SI  material,  and 
that  it  was  compensating  some  impurity  (probably  carbon) 
resulting  from  the  LEC  process. 

Before  considering  the  nature  of  this  compensation  process 
it  is  important  to  ask  a  more  fundamental  question.  Is  this  EL2 
defect  a  process  dependent  effect  that  will  be  eliminated  in  the 
next  few  years  by  better  growth  methods  for  GaAs?  Although  it  is 
likely  that  better  methods  will  be  developed  to  grow  GaAs 
material,  a  vast  body  of  information  exists  for  Czochalski  growth 
for  silicon.  Although  better  growth  techniques  may  be  developed, 
from  an  engineering  and  economy  standpoint  LEC  material  may 
always  be  the  material  of  choice. 
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2.  Origin  of  the  EL2  defect 


In  1984,  Holmes  et.al.  [14]  measured  two-dimensional 
dislocation  density  profiles  and  observed  a  correlation  between 
dislocation  density  and  EL2  concentration.  This  correlation 
had  been  observed  before,  but  only  in  one-dimensional  cross- 
section  [15-17]. 

In  1980,  Jordan  et.al.  [18]  analyzed  two-dimensional 
thermal  stress  profiles  in  a  growing  Czochralski  boule  and 
observed  a  correlation  between  thermal  stress  and  dislocation 
density.  Since  1980,  several  studies  have  verified  this 
correlation  between  theoretical  stress  profiles  and 
experimentally  measured  dislocation  densities  [19-21], 

Clearly,  if  thermal  stress  correlates  to  dislocation 
density,  and  dislocation  density  correlates  to  EL2  concentration, 
then  thermal  stress  correlates  to  EL2  concentration. 

The  KEY  question  is  if  thermal  stress  is  LINKED  to  EL2 
concentration.  This  possibility  was  first  suggested  by  Martin 
et.  al.  in  1981  [22].  Evidence  for  this  link  was  provided  by 
Weber  et.al.  (1982)  when  he  measured  the  creation  of  EL2  defect 
during  plastic  deformation  of  GaAs  [23]. 
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A  variety  of  mechanisms  have  been  proposed  to  account  for 
the  stress  induced  production  of  EL2  in  GaAs .  These  mechanisms 
fall  into  three  principle  categories: 

1.  Stress  forms  dislocations,  and  dislocations  create  EL2 

2.  Stress  forms  EL2,  and  EL2  condenses  into  dislocations 

3.  EL2  and  dislocations  are  independent 

In  1981,  Parsey  et.al.  showed  that  it  is  possible  to  grow 
dislocation-free  material  with  an  EL2  concentration  comparable  to 
conventional  material  [24].  This  result  suggests  that 
dislocations  do  not  create  EL2 . 

In  1981,  Chen  et.al.  showed  that  dislocations  always 
increase  from  seed  to  tail  in  a  CZ-LEC  GaAs  crystal,  but  that  EL2 
concentration  depends  on  the  stoichometry  [25].  This  suggests 
that  EL2  condensing  into  dislocations  is  not  the  principle 
mechanism  for  the  formation  of  dislocations. 

This  evidence  suggests  that  stress  produces  dislocations  and 
the  EL2  defect,  but  that  the  two  processes  are  independent. 
Thus,  prediction  of  the  stress  induced  during  growth  of  the  GaAs 
boule  could  be  used  to  predict  the  EL2  concentrat ion . 


3.  Mathematical  analysis  of  stress  in  a  growing  Czochralski 
boule . 


The  system  being  modeled  is  a  GaAs  boule  during  Czochralski 
growth.  (See  Figure  2)  The  key  item  being  calculated  is  the 
thermal  stress  induced  in  the  boule  as  a  result  of  the  growth 
process.  (Note  that  this  is  NOT  the  residual  stress  contained  in 
the  boule  after  the  boule  is  removed  from  the  melt  and  cooled) 
The  key  assumption  is  that  the  EL2  defect  is  formed  during  the 
growth  process  at  the  time  when  the  thermal  stress  is  maximized. 

Jordan's  1979  paper  [26]  contains  a  careful  analysis  of  heat 
conduction  and  thermal  stress  during  Czochralski  growth.  The  key 
equations  from  this  paper  are  summarized  in  the  following  few 
pages.  (Note:  In  order  to  increase  the  usefulness  of  this 
summary,  the  derivation  is  presented  backwards;  the  answer 
preceeds  the  fundamental  equations  used  to  obtain  it.) 
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Figure  2  -  Czochralski  Growth  of  a  GaAs  Boule 

e  system  being  modeled  is  a  GaAs  boule  during  Czochralski 
owth.  The  key  item  being  calculated  is  the  thermal  stress 
auced  in  the  boule  as  a  result  of  the  growth  process.  (Note 
at  this  is  NOT  the  residual  stress  contained  in  the  boule  after 
e  coule  is  removed  from  the  melt  and  cooled)  The  key  assumption 
that  the  EL2  defect  is  formed  during  the  growth  process  at  the 
me  that  the  thermal  stress  is  maximized. 


Initially,  a  few  boule  specific  parameters  must  be  obtained: 


K  =  THERMAL  CONDUCTIVITY  (WATT/CM  K), 
r  =  EVALUATION  RADIUS  (cm) 

0  =  EVALUATION  ANGLE  (RADIANS) 

P  =  PULL  RATE  (CM/SECOND) 

S  =  Z+PT  =  EVALUATION  LOCATION  (CM) 
a  =  EXPANSION  COEFFICIENT  (1/K) 

T  »  boule  temperature  (K) 

Tr  »  melt  temperature  (K) 

l  =  BORIC  OXIDE  THICKNESS  (CM) 


To  =  BOULE  RADIUS  (cm) 

T  -  PULL  TIME  (SECONDS) 
h  =  HEAT  TRANSFER  (1/CH) 

E  =  Young's  modulus  (dynes/cm2) 
v  =  Poisson's  ratio 

Ta  =  AMBIENT  TEMPERATURE  (K) 
f  =  EMISSIVITY 
P  -  PRESSURE  (DYNES/CM2) 


In  order  to  simplify  analysis,  the  following  additional 
terms  are  defined  from  the  fundamental  information: 


hi  =  hr0 
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Assuming  that  the  stresses  are  mostly  elastic,  the  total 
stress  at  a  given  (rho,  theta)  location  in  a  GaAs  wafer  can  be 


expressed  as  the  sum  of  the  stresses  along  all  slip  planes. 
<?iot  =  4|<t/|  +  2|<t//|  +  2|<r///|  -f  2|<t/v|  +  2|<tv| 


In  the  GaAs  system  (analyzed  quasi-isotropically) ,  there  are 
twelve  possible  slip  planes.  Five  of  these  slip  planes  are 
independent.  The  stresses  on  these  five  independent  slip  planes 
can  be  expressed  as  follows: 


direction/plane 

[-110]  /  (111) 
[110]  /  (-111) 
[1-10]  /  (-1-11) 
[-1-10]  /  (l-ll) 


[0-11]  /  (ill) 

[Oil]  /  (-1-11) 


[0-11]  /  (-111) 
[Oil]  /  (1-11) 


[10-1]  /  (111) 
[-10-1]  /  (-1-11) 


[-10-1]  /  (1-11) 

[10-1]  /  (1-11) 
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Where,  in  these  equations  the  following  alternative  radial 


and  axial  stress  terms  have  been  used: 


Or  =  Or  ~ 
cTz  —  oz-  OQ 


The  radial,  angular  and  axial  stress  terms  contained  in  the 
alternative  stress  terms  can  be  obtained  from: 


V  (  Vi»Mn  (*e  -  *)  +  ftn  COShjffn  (V  t  ~  *) 

“  \(/l5  +  «?J  J0  (««)  (^Psinl‘At^<  +  Pn  cash  t) 


j  1  (an)  h  {(*np) 
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Where  See  is  given  by: 


See  =  2/1-——  (Tf  -  Ta)exp*'*'2r0 

l-i/ 

There  are  two  possible  solutions  for  the  thermal  stress  in  the 


growing  boule. 


The  previous  solution  results  from  the 


assumption  of  zero  axial  strain  on  the  growing  boule.  However, 
if  zero  axial  force  is  assumed,  then  the  final  solution  differs 
in  the  axial  stress  term.  The  axial  stress  term  for  zero  axial 
force  is  given  by: 


<JT  =  Sr 
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It  is  still  unclear  which  of  the  two  assumptions  (zero  axial 
strain,  or  zero  axial  force)  is  the  better  model  of  the  growing 
Czochralski  boule. 


All  of  the  terms  in  the  equations  on  the  previous  page  are 
known  material  parameters  EXCEPT  for  the  terms  alpha  and  beta. 
Beta  is  a  function  of  alpha,  defined  by: 


0n  = 


1/2 


However,  alpha  must  be  obtained  from  the  characteristic 
equation: 

-anJ\  (On)  +  hiJo  (on)  =  0 


The  radial,  angular  and  axial  stress  terms  detailed 
previously  are  derived  by  substituting  the  temperature  profile 
for  the  growing  boule  into  the  fundamental  integral  stress 
equations  for  cylindrical  geometry.  The  fundamental  integral 
stress  equations  (assuming  zero  axial  strain)  are: 


(Tr  = 


Trdr- 
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With  the  assumption  of  zero  axial  force,  the  integral 
equation  for  the  axial  stress  becomes: 


The  temperature  profile  required  in  the  integral  stress 
equations  is  derived  from  the  fundamental  heat  conduction 
equation  for  a  cylindrical  geometry. 

ar  (d2T  l  dT  aJn 

3t  K\dr2  +  r  dr  dzx) 

This  fundamental  heat  conduction  equation  can  be  simplified 
into  the  quasi-steady  state  form,  and  written  with  the 
dimensionless  terms  rho,  eta  and  psi: 

d2T  \  dT  d2T  dT 

- r  + - 1 - 7  =  p  i - 

dp  p  dp  341  d^p 

The  following  product  function  can  be  used  to  separate  the 
variables : 

T  =  ep'+/2fUp)  +  {4,) 
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So,  the  final  solution  to  the  cylindrical  quasi-steady  state 


thermal  conduction  equation  is: 


T  =  Ta  +  ep''¥n  X  Jo(a„p)  ~  (Vipsinh  /?„(«//,  -  if,) 


B„ 


ft* 


+  /ff*cosh  //..(ft  -  v//)], 


The  boundary  conditions  for  the  growing  boule  are: 


+  hi(T  0 

ay 


dT 


—  +  h(T-T..)\,-P<  =  Q 
ds 


These  can  be  applied  to  the  solution  of  the  cylindrical 
quasi-steady  state  thermal  conduction  equation,  to  yield: 


dJu 


dp 


( cr„p )  +  h\ Jo{a.,p)  |„-i  =  0 


Applying  the  Bessel  function  recursion  relation: 


to  the  previous  equation,  results  in  the  Bessel  funrt  ion 

characteristic  equation  which  we  showed  previously. 


+  /iiJo(Qfo)  =  0 


4.  Behaviour  of  the  Stress  in  a  Growing  CZ-LEC  Boule 


Several  of  the  growth  parameters  are  key  to  understanding 
the  behaviour  of  the  growing  Czochralski  boule.  These  are: 

h  =  The  heat  transfer  coefficient  between  the  boule  and 
the  ambient  environment 

ro  =  The  radius  of  the  boule 

p  =  The  growth  rate  of  the  boule 

t  =  The  time  to  grow  the  boule 

s  =  z  +  pt  «=  the  distance  from  the  melt  interface  to  the 
location  where  the  stress  is  being  calculated. 

"h"  governs  the  heat  transfer  between  the  surface  of  the 
boule  and  the  ambient  environment.  (Less  heat  transfer 
smaller  "h")  Intuitively,  we  would  expect  reduced  stress  for 
reduced  heat  transfer,  so  reducing  "h"  should  reduce  the  stress. 
This  intuitive  conclusion  is  supported  by  Figure  3,  where  the 
normalized  stress  is  compared  for  h  =  1.8  and  h  =  0.6.  (The 
zero  axial  strain  solution  is  assumed,  and  other  parameters  are 
held  constant  at:  p=1.94e-4  cm/sec,  t=51546  sec,  r=3.8i  cm, 

s=l . 5  cm  and  theta  =  0.785  radians) 

"ro"  is  the  radius  of  the  final  boule.  Intuitively,  v.'o 
v/ould  expect  that  the  smaller  the  boule,  the  greater  the  stress 
involved  during  growth.  This  is  borne  out  by  Figure  4,  vh ■  r. 


I, 
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compares  the  normalized  stress  at  theta  =  .785  for  ro  =  2.54  cm 
(2"  diameter  boule)  and  ro  =  3.81  (3”  diameter  boule) .  (The 
zero  axial  strain  solution  is  assumed,  and  other  parameters  are 
held  constant  at:  p=1.94e-4  cm/sec,  t=51546  sec,  h  =  0.6, 
s=l . 5  cm  and  theta  =  0.785  radians) 


d 
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"p"  is  the  pull  rate  and  "t"  is  the  amount  of  time  it  takes 
to  pull  the  boule.  Intuitively,  we  would  expect  that  the  longer 
it  takes  to  pull  the  boule  the  less  the  stress  will  be.  This  is 
true,  as  Figure  5  shows,  but  the  effect  is  not  as  significant  as 
might  be  expected.  (The  zero  axial  strain  solution  is  assumed, 
and  the  other  parameters  are  held  constant  at  a  constant  length 
of  10  cm:  h  *  0.6,  ro  «  3.81  cm,  s=l.o  cm  and  theta  =  0.785 
radians) 

"s"  is  the  distance  from  the  melt/boule  interface  to  the 
location  where  the  stress  will  be  computed.  Clearly,  as  this 
distance  increases,  the  stress  will  drop.  This  intuitive 
conclusion  is  borne  out  by  the  results  in  Figure  6,  where  the 
normalized  stress  is  compared  as  a  function  of  "s".  (The  zero 
axial  strain  solution  is  assumed,  and  the  other  parameters  are 
held  constant  at:  h  =  0.6,  ro  =  3.81  cm,  p=1.94e-4,  t  - '  1  v;  h  ; 
theta  =  0.785  radians) 
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Mathematically ,  the  only  difference  between  the  zero  axial 
strain  and  the  zero  axial  force  boundary  conditions  is  a  factor 
of  Poisson's  ratio  in  the  axial  stress  equation.  However,  the 
total  stress  results  differ  significantly  depending  on  the  chosen 
assumption.  Total  stress  results  are  shown  in  Figure  7  for  both 
the  zero  axial  strain  and  the  zero  axial  force  assumptions. 
Note  that  the  assumption  of  zero  axial  force  results  in  greater 
overall  stresses,  as  well  as  a  stress  minimum  placed  at  a  greater 
distance  from  the  wafer  center  than  tlii  zero  axial  strain 
assumption. 


Normalized  Stress 


EFFECT  OF  BOULE  DIAMETER 
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Figure  4  -  Effect  of  Boule  Diameter 

"ro"  is  the  radius  of  the  final  boule.  Intuitively,  we  would 
expect  that  the  smaller  the  boule,  the  greater  the  stress 
involved  during  growth.  This  is  borne  out  by  the  above  graph 
v/nich  compares  the  normalized  stress  at  theta  =  .785  for  ro 
:  .  r  cm  (2"  diameter  boule)  and  ro  =  3.81  (3"  diameter  bou.c)  . 

:r.o  zero  axial  strain  solution  is  assumed,  and  other  paramete 
are  held  constant  at:  p-l.g4o-4  cm/sec,  t=51546  sec,  h  =  0. 
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Figure  5  -  Effect  of  Pull  Rate  and  Time 

"p"  is  the  pull  rate  and  "t"  is  the  amount  of  time  it  ta.:e 
pull  the  boule.  Intuitively,  we  would  expect  that  the  long 
takes  to  pull  the  boule  the  less  the  stress  will  be.  Thi 
true,  as  Figure  5  shows,  but  the  effect  is  not  as  signifies 
might  be  expected.  (The  zero  axial  strain  solution  is  ass 
and  the  other  parameters  are  held  constant  at  a  constant  1 
of  10  cm:  h  =  0.6,  ro  =  3.81  cm,  s=1.0  cm  and  theta  = 

radians ) 
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Figure  7  -  Zero  Axial  Force  versus  Zero  Axial  Strain 

Mathematically,  the  only  difference  between  the  zero  axial  strain 
and  the  zero  axial  force  boundary  conditions  is  a  factor  of 
Poisson's  ratio  in  the  axial  stress  equation.  However,  the 
total  stress  results  differ  significantly  depending  on  the  chosen 
assumption.  Total  stress  results  are  shown  above  for  both  the 
zero  axial  strain  and  the  zero  axial  force  assumptions.  Note 
that  the  assumption  of  zero  axial  force  results  in  greater 
overall  stresses,  as  well  as  a  stress  minimum  placed  at  a  greater 
distance  from  the  wafer  center  than  the  zero  axial  strain 
assumption.  (The  other  parameters  are  held  constant  at:  h  =  0.6, 
ro  =  3.81  cm,  p=l.94e-4  cm/sec,  t=51546  sec,  s  =  1  cm  and  theta  = 
0.785  radians) 


5. 


Obtaining  the  stress  conversion  coefficient 


Jordan's  equations  are  usually  solved  for  a  "normalized- 
stress"  which  differs  from  the  actual  stress  by  a  factor  See  (see 
previous) .  Ideally,  the  EL2  concentration  should  be  obtained  by 
multiplying  this  normalized  stress  by  a  factor  Mcc  which  consists 
of  the  product  of  See  and  the  stress-conversion-coefficient. 
This  factor  Mcc  is  constant  for  all  boules  grown  with  identical 
temperature  distributions  and  boric  oxide  thicknesses,  (the 
radii  and  pull  rates  may  differ) .  Note  that  this  method  of  EL2 
calculation  does  not  require  knowledge  of  the  values  of  T,  Ta  and 
Tf.  This  is  advantageous  since  these  temperatures  are  difficult 
to  measure  experimentally. 

We  empirically  determined  Mcc  by  using  a  computer  model 
(developed  from  Jordan's  stress  equations)  to  compare  theoretical 
stress  values  (both  for  the  assumption  of  zero  axial  strain  and 
of  zero  axial  force)  with  experimental  results  from  Brozel  et.al. 
[27].  The  results  of  this  comparison  are  shown  in  Figure  8. 
The  zero  axial  strain  assumption  yields  a  better  match  between 
theory  and  experiment  with  an  empirically  determined  value  for 
Mcc  obtained  from  this  comparison  is  9el6  [EL2]/cm4. 
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Figure  8  -  Obtaining  the  Stress  Conversion  Coefficient 

We  empirically  determined  Mcc  by  using  a  computer  model 
(developed  from  Jordan’s  stress  equations)  to  compare  theoretical 
stress  values  (both  for  the  assumption  of  zero  axial  strain  and 
of  zero  axial  force)  with  experimental  results  from  Brozel  et.al. 
[27].  The  results  of  this  comparison  are  shown  in  Figure  8. 
The  zero  axial  strain  assumption  yields  a  better  match  between 
theory  and  experiment  with  an  empirically  determined  value  for 
Mcc  obtained  from  this  comparison  is  9el6  [EL2]/cm4.  (The  other 
parameters  are  held  constant  at:  h  =  0.6,  ro  =  3.81  cm,  p=i.94e-4 
cm/sec,  t=5 154 6  sec,  s  =  1  cm  and  theta  =  0.785  radians) 


v.  • 


6.  Predicting  the  two-dimensional  EL2  distribution 

Once  the  "stress-conversion-coefficient"  is  known,  it  is 
possible  to  predict  two-dimensional  EL2  concentration  profiles 
using  the  previously  described  equations.  We  developed  a 
computer  model  capable  of  calculating  the  EL2  distributions  from 
fundamental  material  parameters.  Figure  9  shows  a  comparison 
between  the  predicted  EL2  profiles  (for  both  the  zero  axial 
strain  and  the  zero  axial  force  assumptions)  and  experimental 
profiles  (measured  by  Holmes) .  Figures  10  and  11  show  detailed 
views  of  the  predicted  EL2  profiles  for  each  axial  assumption. 
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Figure  9  -  Two  Dimensional  EL2  Concentration  Profiles 

Once  the  "stress-conversion-coefficient"  is  known,  it  is  possible 
to  predict  two-dimensional  EL2  concentration  profiles  using  the 
previously  described  equations.  We  developed  a  computer  model 
capable  of  calculating  the  EL2  distributions  from  fundamental 
material  parameters.  The  above  figure  shows  a  comparison 
between  the  predicted  EL2  profiles  (using  both  the  zero  axial 
stress  and  zero  axial  strain  assumptions)  and  experimental 
profiles  (measured  by  Holmes) .  (The  other  parameters  are  held 
constant  at:  h  =  0.6,  ro  =  3.81  cm,  p=1.94e-4  cm/sec,  t=51546  sec, 
s  =  1  cm  and  theta  =  0.785  radians) 
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EL2  CONCENTRATION  IN  10**15  PER  CC 
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Figure  10  -  The  Zero  Axial  Strain  Assumption 

The  above  figure  shows  a  detailed  view  of  the  EL2  conce: 
profile  resulting  from  using  the  zero  axial  strain  ass; 
(The  other  parameters  are  held  constant  at:  h  =  0.6,  ro  = 
p=1.94e-4  cm/sec,  t=51546  sec,  s  =  0.6  cm  and  theta 
radians) 
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Figure  11  -  The  Zero  Axial  Force  Assumption 

The  above  figure  shows  a  detailed  view  of  the  FL2  concentratic 
profile  resulting  from  using  the  zero  axial  force  assumption 
CThe  other  parameters  are  held  constant  at:  h  =  0.6,  ro =  3.61  cr 
p=1.94e-4  cm/sec,  t=51546  sec,  s  =  0.6  cm  and  theta  =  0.7P 
radians ) 
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IV.  ANNEALING  THE  EL2  DEFECT 

1.  Introduction 

The  critical  step  in  fabricating  an  ion-implanted  MESFET 
device  is  the  activation  anneal.  It  is  vitally  important  that 
the  electrical  characteristics  of  the  GaAs  substrate  not  change 
during  the  activation  anneal. 

We  developed  a  quasi-two-dimensional  computer  model  to 
simulate  the  effects  of  an  activation  anneal  on  the  EL2  defect 
profile.  This  model  takes  the  theoretical  EL2  defect  density  and 
performs  a  one-dimensional  Gaussian  diffusion  at  each  calculated 
point  in  the  substrate. 

Initially,  the  approximate  diffusion  coefficient  for  the  EL2 
defect  in  GaAs  was  determined  by  taking  existing  experimental 
results  and  curve-fitting  them  to  the  results  of  the  computer 
model.  Once  the  diffusion  coefficient  was  obtained,  it  was  used 
in  conjunction  with  the  computer  model  to  predict  EL2  defect 
density  profiles  before  and  after  annealing. 
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2.  Initial  Observation  of  Annealing 

A  comparison  between  theory  and  experiment  was  performed 
using  data  from  Holmes  [28].  The  results  of  this  comparison  tor 
both  the  zero  axial  strain  and  zero  axial  force  assumptions 
(using  the  multiplicative  coefficient  of  1.9  el7  [EL2]/cmA4)  are 
shown  in  Figure  12.  Although  the  overall  agreement  between 
theory  and  experiment  is  quite  good  for  the  zero  axial  strain 
assumption,  the  edges  do  deviate  considerably  from  theory. 

We  postulated  that  this  deviation  was  due  to  a  post-growth 
annealing  process  occurring  as  a  result  of  the  boule  cooling  in 
the  puller  for  many  hours.  We  therefore  expanded  the  computer 
model  to  take  theoretical  EL2  concentrations  (as  predicted  by  the 
stress  calculations)  and  diffuse  these  throughout  the  boule. 


Figure  12  -  Comparison  of  Theory  with  Experiment  (Holmes) 

A  comparison  between  theory  and  experiment  was  performed  us 
data  from  Holmes  [28].  The  zero  axial  strain  assumption  yield 
good  correlation  between  theory  and  experiment  with 

empirically  determined  value  for  Mcc  of  1.9el5  [EL2],'cm4.  < 

other  parameters  are  held  constant  at:  h  ^  0.6,  ro  =  2.01 

p=1.94e-4  cm/sec,  t=51546  sec,  s  =1  cm  and  theta  =  o . 
radians) 


Determination  of  the  Approximate  Diffusion  Coefficient 


3  . 

Figure  9  shows  experimental  annealing  results  from  Holme;- 
[24].  Using  our  annealing  model,  we  roughly  evaluated  the 
diffusion  constant  of  EL2  at  950  C  by  comparing  theoretically 
calculated  annealed  profiles  against  the  actual  anneale.i 
profiles.  These  calculations  result  in  a  diffusion  constan- 

of  7  (+-  5)  e-7  cm2/sec. 

The  value  obtained  for  this  diffusion  constant  is  quite 
high.  The  magnitude  of  this  value  suggests  that  EL2  is  diffusing 


by  an  interstitial 

mechanism  as  opposed 

to  a 

subst it iona 1 

mechanism. 

This,  in 

turn,  suggests  that 

some 

part 

of  the 

mechanism 

responsible 

for  generating  EL2  may  be 

based 

on  an 

interstitial  defect. 


These  annealing 

results 

suggest 

the  following  model 

During  Czochralski  growth,  thcr 

rr.al  at  re 

sr.es  cause  modi  f  ioat 

~£  the  electric  field 

in  the  ionic  GaAs 

crystal  . 
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modifications  create 

charged 
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attract ive 
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interstitial  arsenic 
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anneal ing , 

the  ex.-er 

interstitial  arsenics  become  more  mobile,  and  migrate  to  oth‘o- 
charged  sites  in  the  crystal. 
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Figure  13  -  Annealing  the  EL2  Defect  (Holmes) 

This  graph  shows  experimental  results  from  Holmes  obtained  from 
annealing  the  EL2  defect  in  GaAs  [24].  Using  our  annealing 
model,  we  roughly  evaluated  the  diffusion  constant  of  EL2  at  950 
C  by  comparing  theoretically  calculated  annealed  profiles  against 
he  actual  annealed  profiles.  These  calculations  result  in  a 

iffusion  constant  of  7  (+-  5)  e-7  cm2/sec. 


4.  Comparison  of  the  Annealing  Model  with  Experiment 


As  we  mentioned  earlier,  we  postulated  that  the  difference 
in  curve  shape  between  the  theoretical  results  and  Holmes's 
experimental  results  was  due  to  a  post-growth  annealing  process 
occurring  as  a  result  of  the  boule  cooling  in  the  puller  for 
many  hours.  Using  this  assumption,  as  well  as  assuming 
outdif fusion  from  the  GaAs  boule,  we  recalculated  the  EL2 
distribution. 

Figures  14  compares  the  diffused  EL2  results  with  the 
experimental  values.  The  correlation  between  the  diffused  EL2 
results  and  the  actual  experimental  values  (from  Holmes)  is  quite 
good.  This  suggests  that  post-growth  whole-boule  anneal  may 
alter  the  EL2  concentration  results  as  obtained  from  stress 
calculations.  It  also  suggests  that  post-growth  whole-boule 
anneal  may  be  a  way  to  equalize  EL2  concentration  across  the 
boule . 
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ANNEALING  THE  EL2  DEFECT 

(Experimental  Results  from  Holmes) 
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Figure  M  -  Comparison  of  annealed  and  unannealed  results  for 

the  EL2  defect 

This  graph  compares  the  diffused  EL2  results  (obtained  from  the 
annealing  model  as  applied  to  the  stress  analysis)  with  the 
experimental  values  from  Holmes. 
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ELECTRICAL  CHARACTERISTICS  OF  THE  SI-GaAs  SUBSTRATE 
AND  THE  THERMAL  STABILITY  OF  THESE  CHARACTERISTICS 


1.  Introduction 


Clearly,  the  variation  of  the  EL2  defect  across  the 
substrate  is  unimportant  if  it  does  not  affect  device 
performance.  Since  substrate  resistivity  is  of  primary  concern 
in  an  ion-implanted  device,  our  first  objective  was  to  predict 
the  resistivity  at  each  location  on  the  wafer  as  a  function  of 
the  EL2  concentration. 


A  Fermi-level  computer  model  was  constructed  to  calculate 
the  Fermi-level  at  each  point  in  the  wafer  as  a  function  of  four 
arbitrary  impurities.  Once  obtained,  the  Fermi  level  was  then 
used  to  calculate  the  resistivity  at  each  location  in  the  wafer. 
Theoretical  results  were  obtained  for  unannealed  material 
containing  EL2 ,  carbon,  iron  and  silicon. 
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2.  Theory  of  the  Fermi-level  Analysis 


There  are  a  variety  of  impurities  and  defects  in  SI-LEC 


material  (see  Figure  15) 


The  most  important  are  carbon 


(acceptor) ,  iron  (acceptor) ,  EL2  (donor)  and  silicon  (donor) . 

The  relative  concentrations  and  energies  of  these  various 
impurities  determines  the  location  of  the  Fermi-level  in  the 
bandgap.  Once  calculated,  the  location  of  the  Fermi-level  in 
the  gap  can  be  used  to  obtain  many  other  electrical  properties  of 
interest  (resistivity,  and  trap-filled  limit  voltage  for 
example) . 


The  charge  neutrality  condition  for  a  substrate  containing 
four  major  impurity  levels  is: 


,  v  (~)  (“)  (  +  )  (  +  ) 

71  ■  1  A-tron  +  PiA-carbon  =  P  +  ND_£L2  +  ND-silit 


The  electron  and  hole  density  (n  and  p)  are  given  in  terms 
of  the  Fermi-level  as: 
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Figure  15  -  Summary  of  Defects  and  Impurities  in  GaAs 

There  are  a  variety  of  impurities  and  defects  in  SI-LEC  material, 
the  most  important  are  carbon  (acceptor) ,  iron  (acceptor) ,  EL2 
(donor)  and  silicon  (donor) .  The  relative  concentrations  and 
energies  of  these  various  impurities  determines  the  location  of 
the  Fermi-level  in  the  bandgap.  Once  calculated,  the  location 
of  the  Fermi-level  in  the  gap  can  be  used  to  obtain  many  other 
electrical  properties  of  interest  (resistivity,  and  trap-filled 
limit  voltage  for  example) . 


The  concentrations  of  the  various  defects  and  impurities  are 
also  given  in  terms  of  the  Fermi-level  as: 


Na  = 


1+8  ex'»(SakT&) 


No  =  No  1- 


,+?exp( 
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The  concentrations  of  carbon,  iron  and  silicon  are 
essentially  constant  across  an  SI-LEC  wafer.  Unfortunately,  the 


EL2  concentration  displays  a  characteristic  W  pattern. 


This 


means  that  the  Fermi-level  must  be  recalculated  for  each  point  on 
the  wafer. 


Once  the  Fermi-level  is  determined  for  each  point,  then  it 
can  be  used  to  calculate  various  electrical  parameters  of 
interest  at  that  point.  For  example,  the  electron  and  hole 
densities  (n  and  p)  can  be  calculated  from  the  Fermi-level: 

n  =  N^p(—CkTEf)) 
p  =  N„exp  (~( E,k~E 

These  can  then  be  used  to  determine  the  resistivity  on  a 
point  by  point  basis. 
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3.  Calculation  for  a  Wafer 

Holmes's  data  [30]  was  selected  for  an  initial  computation 
of  the  Fermi-level  and  resistivity  in  a  real  system.  The 
characteristics  of  the  boule  and  wafer  are: 

1.  Boule  pull  time  =  51546  seconds 

2.  Boule  pull  rate  =  1.94  e-4  cm/sec 

3.  Assumed  Mcc  =*  1.7  el7  [EL2]/cm4 

4.  Carbon  *  6  el5  [C]/cm3 

5.  Iron  *  3  el5  [Fe]/cm3 

6.  Silicon  -  2  el5  [Si]/cm3 

First,  the  computer  model  was  run  with  these  values  to 
create  a  database  containing  the  EL2  concentration  as  a  function 
of  location.  Then,  the  Fermi-level  and  resistivity 

calculated.  These  results  are  shown  in  Figures  16-17. 


were 


Figure  16  -  Computation  of  the  Fermi-level 

This  graph  shows  an  initial  computation  of  the  Fermi-level  in  a 
real  system.  The  system  parameters  are:  Boule  pull  time  =  51546 
seconds,  Boule  pull  rate  ■  1.94  e-4  cm/sec,  h  *  0.6,  s  =  1  cm,  r 
=  3.81  cm,  assumed  Mcc  -  1.7  el7  [EL2]/cm4,  Carbon  =  6  el5 
[Cj/cm3,  Iron  =•  3  el5  [Fe]/cm3  and  Silicon  -  2el5  [Si]/cm3. 
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Figure  17  -  Computation  of  the  Resistivity 

This  graph  shows  an  initial  computation  of  the  resistivity  in  a 
real  system.  The  system  parameters  are:  Boule  pull  time  =  51546 
seconds,  Boule  pull  rate  =  1.94  e-4  cm/sec,  h  »  0.6,  s  =  1  cm,  r 
“  3.81  cm,  assumed  Mcc  =  1.7  el7  [EL2]/cm4,  Carbon  =  6  el5 

[C]/cm3,  Iron  =  3  el5  [Fe]/cm3  and  Silicon  =  2  el5  [Si]/cm3. 


4 .  Summary 


The  previous  results  are  extremely  interesting.  As  we 
might  expect,  the  extremes  of  the  EL2  concentration  curve  display 
the  greatest  deviation  in  Fermi-level  and  resistivity.  What 
this  means  is  that  an  improvement  in  the  uniformity  of  the  EL2 
distribution  results  in  an  improvement  in  the  uniformity  of  the 
substrate  resistivity.  Since  annealing  improves  the  uniformity 
of  the  EL2  distribution,  it  is  clear  that  the  quality  of  the  GaAs 
substrate  can  be  improved  by  whole-boule  post-growth  annealing. 


VI. CONCLUSIONS 


The  key  technical  issue  in  the  growth  of  SI-LEC  material  is 


understanding  and  controlling  the  EL2  anti-site  defect. 


With 


the  results  from  this  project  in  mind,  answers  can  be  proposed  to 
the  six  questions  introduced  earlier  in  this  report: 

1.  What  is  the  physical  mechanism  which  creates  EL2? 

The  evidence  strongly  suggests  that  thermal  stress 
occurring  during  growth  of  the  Czochralski  boule  creates 
EL2.  EL2  distributions  can  be  accurately  predicted  by 
calculating  the  theoretical  stress  and  applying  a  simple 
multiplicative  "stress-conversion-coefficient" . 

2.  Can  the  EL2  concentration  profile  be  predicted  from 

fundamental  information? 

The  EL2  concentration  profiles  can  be  predicted  by 
applying  a  "stress-conversion-coefficient"  to  the 
results  of  a  theoretical  thermal  stress  analysis.  The 
theoretical  analysis  uses  only  fundamental  material  and 
growth  parameters,  but  the  "stress-conversion- 
coefficient"  is  determined  empirically.  Ideally,  it 
would  be  desirable  to  also  calculate  the  "stress- 
conversion-coefficient"  from  fundamental  material  and 
growth  parameters. 

3.  What  is  the  origin  of  the  ubiquitious  "W"  shape  for 
radial  EL2  density  profiles? 

GaAs  possesses  a  zincblende  structure.  Under  stress, 
materials  with  this  crystal  structure  display  a  two- 
dimensional  stress  pattern  which  resembles  a  four-loaf 
clover.  If  the  stress  is  plotted  along  a  diameter,  it 
will  possess  a  characteristic  "W"  shape.  Since  the 
stress  is  directly  related  to  the  EL2  concentration, 
then  the  EL2  profile  is  also  "W"  shaped. 
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4.  Can  the  MW"-shaped  EL2  concentration  profile  be  made 
constant? 

EL2  diffusing  in  GaAs  possesses  a  rather  high  diffusion 
coefficient.  Thus,  whole  boule  annealing  of  GaAs  can 
be  used  to  reduce  the  variation  in  EL2  concentration. 


5.  Does  the  compensation  between  the  "W"  shaped  EL2 
concentration  profile  and  the  constant  carbon  concentration 
result  in  variation  in  electrical  parameters? 

As  one  might  expect,  the  extremes  of  the  EL2 
concentration  curve  display  the  greatest  deviation  in 
Fermi-level  and  resistivity. 


6.  Is  the  compensation  between  EL2  and  carbon  thermally  stable? 

The  only  effect  of  long  term  (greater  than  one  hour) 
thermal  anneal  is  to  reduce  the  variations  in  the  EL2 
density  profile.  Thus,  not  only  is  the  compensation 
thermally  stable,  but  the  material  tends  to  get  better 
as  it  is  annealed. 


Semi-insulating  Gallium  Arsenide  substrates  offer  the 
advantages  of  device  isolation,  low  interconnect  capacitance, 
simplified  processing,  and  direct  ion-implantation  for  device 
fabrication.  Hopefully  the  results  of  this  project  will  aid  in 
improving  the  quality  and  homogeneity  of  existing  SI-LEC  material 
in  order  to  help  achieve  the  tremendous  advantages  of  GaAs 
technology. 
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